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RELATIVE MOTION IN THE TERMINAL PHASE OF INTERCEPTION CEP 
A SATELLITE OR A BALLISTIC MISSILE 
By Richard A. Hord 

SUMMARY 


The interception of a satellite or a ballistic missile essentially 
outside the earth's atmosphere is considered with emphasis on the ter- 
minal phase of the attack vehicle's flight. For the discussion of the 
relative motion, a coordinate system is selected which has its origin 
at the target and axes always parallel to lines fixed in an inertial 
frame. The resulting vector equation of motion contains, in addition 
to the thrust term, an apparent gravitational term. Bounds on the magni- 
tude of this apparent acceleration of the interceptor are obtained which 
permit gaging its importance in relation to thrust acceleration in the 
terminal phase of interception. Since the apparent gravitational accel- 
eration approaches zero as the interceptor approaches the target, pre- 
liminary analyses of the terminal phase of interception can generally be 
carried out by neglecting gravity altogether. Two types of interception 
are discussed mathematically and a numerical example is worked out for 
each type . The type of interception for which the terminal relative 
approach speed is small compared with the target's speed is, within the 
limits of the present study, found to be feasible. The practicability 
of head-on interception with mass dispersal prior to impact is more 
difficult to assess; in the example considered, the mass required to be 
dispersed is found to be of the order of magnitude of the mass of the 
target. 


INTRODUCTION 


The use of pilotless aircraft in the destructive interception of 
an airplane is presently a well-developed field of engineering both in 
theory and in practice. Ihe corresponding literature is extensive and 
no attempt need be made to review it here. (See, for example, ref. 1.) 

The interception of satellites and ballistic missiles, on the other 
hand, remains practically unexplored. This case is characterized by 
speeds so high as to approach escape speed, propulsion and control by 
rockets, and a shift of emphasis, particularly in the terminal phase of 
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flight, from aerodynamic forces to gravitational forces. In addition, 
the motion of a satellite or a "ballistic missile can "be accurately pre- 
dicted over large distances from knowledge of its earlier motion. 

In this paper attention is centered on the terminal phase of the 
interceptor's flight. Equally important problems of initial detection, 
tracking, launching, initial guidance and propulsion, and stabilization, 
for example, are not considered to an appreciable extent in the present 
study. 


SYMBOLS 


a* position vector of attack vehicle in a coordinate system 

with origin at target's center of mass and axes always 
parallel to lines fixed in an inertial frame 

target's frontal area 

b radius of sphere over which dispersed mass is distributed 

c navigational correction coefficient 

— > 

F thrust vector 



m 

m' 

M 

m s 

“md 


vector sum of perturbing forces at r 

vector sum of perturbing forces at p* 

Newton' s universal gravitational constant, 

6.670 x 10" 11 newton-meter 2 /kg 2 = 3.438 x 10 "® lb-ft^slug 2 

mass of attack vehicle 

mass of target 

mass of earth 

mass which strikes target 

minimum destructive mass 


position vector of attack vehicle in a coordinate system 
with origin at earth's center and axes always parallel 
to lines fixed in an inertial frame 


r 
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t time 

U potential function for tide-generating acceleration, or for 

— » 

r 

v^ magnitude of velocity vector da^dt at time t^ 

v ex exit speed of propulsive exhaust gases 

w^ radial speed of dispersed mass 


-» 

r 


ap 


aT 


pr 


9 


apparent gravitational acceleration acting on attack vehicle 

maximum error (maximum distance of center of mass of dispersed 
matter from target's flight path at interception) 

angle which velocity vector da*/dt makes with negative of 
position vector a* at time t’^ 

angle between positive directions of vectors 

angle between positive directions of vectors 

angle between positive directions of vectors 

position vector of target in a coordinate system with origin 
at earth's center and axes always parallel to lines fixed 
in an inertial frame 

radius of earth, 6,371 Inn = 3 >959 miles 
dimensionless correction constant, c^mv^ 
longitude in a spherical coordinate system 


a* and p* 
— » 

a and r 

— » — » 

p and r 


Subscripts : 


1 refers to start of terminal phase of interception 

i refers to start of mass dispersal (type II) 


2 


refers to time at which target is hit 
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GENERAL CONSIDERATIONS 


The problem, of intercepting a satellite vehicle or a ballistic mis- 
sile essentially outside the earth's atmosphere has two features which 
alone are sufficient to set it apart from the interception problems 
encountered in utilizing a guided missile to destroy a conventional air- 
craft. These features are: 

(1) The predictability of the target's flight path over large dis- 
tances. The accurate measurement of the position and velocity of the 
target at some point in time makes possible the computation of its posi- 
tion and velocity at future times by the methods of celestial mechanics. 
This means that the interception program can.be primarily predetermined, 
the subsequent corrections required in flight being comparatively small. 
For that portion of the attack vehicle 1 s flight path which lies in the 
atmosphere, aerodynamic controls can be employed, whereas outside the 
atmosphere propulsive reaction devices must be used for control. 

(2) Hie target's extremely high speed. Because of this, the choice 
of interception program is more critically dependent upon the capabilities 
of the components of the interception system concerned with initial detec- 
tion, tracking, guidance, propulsion, stabilization, and destruction. In 
particular, for an interception system requiring a direct hit, the termi- 
nal approach speed, relative to the target, must be low enough to permit 
the system to overcome its limitations in attack-vehicle maneuverability 
and in component accuracy and response time. An important case is the 
extreme case in which the terminal relative approach speed is small, 

that is, when the attack vehicle's velocity vector closely approximates 
the target's velocity vector immediately prior to the hit. Interception 
programs with this characteristic will be called type I. On the other 
hand, the target's extremely high speed can clearly be used to advantage 
In destroying it by a head-on attack. In this case the demands placed 
upon the interception system probably preclude scoring a direct hit. 
However, the timed dispersal of all or a part of the attack vehicle 1 s 
mass prior to interception can be utilized to reduce the accuracy require- 
ment. A simple and important case is that in which the center of mass of 
the attack vehicle, during and somewhat before dispersal, is in the pro- 
jected orbit of the target. At contact, then, the velocities are equal 
in magnitude and oppositely directed. Head-on interception programs of 
this nature will be called type II. 

The attack vehicle can be launched from the ground or from another 
vehicle such as an aircraft, a satellite, or. .a satelloid. Flight paths 
of ground- launched interceptors of the two types described are shown 
schematically in figure 1. The chief differences between the two flight 
paths are the following: 
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(1) The terminal portion of type I will generally he extended because 
of the small relative approach speed. Thus, the time in flight tends to 
be greater than in the case of type II. 

(2) The launching point of type I will, in general, lie nearer to the 
point of initial detection. This suggests, for example, a single location 
for the equipment for detection, tracking, launching, and so forth, in the 
case of type I. 

(3) Because of the difference in flight-path orientation, the earth's 
rotational motion may be usable to better advantage in one case than in 
the other. 

(4) If considerable variation in point of destruction is allowable, 
the time at which the interceptor is launched need not, in general, be 
as precise for type II as for type I. 

(5) If the eccentricity of the target's elliptical orbit is large 
and if the point of destruction is near the minor axis of the ellipse, 
then, depending upon which way the target’s velocity vector is inclined 
to the horizontal, the time in flight tends to differ considerably for 
the two types. 

These differences, of .course, become less important as the overall 
performance of the interception system increases. 


EQUATIONS OF MOTION 


The equations of motion for the two types are basically the same 
and, outside the earth's atmosphere, may be expressed in vector form as 



-0 ^ ?+ i? 

3 m 


+ - ft. 

m P 


( 1 ) 


— > I 

where r is the position vector of the center of mass of the attack 
vehicle at time t, G is the -universal gravitational constant, M is 

the earth's mass, m is the instantaneous mass (nonincreasing with time), 
—■> — > 

F is the thrust, and F is the vector sum of perturbing forces such as 

P 

those due to the earth's oblateness and motion around the sun, the moon's 

■'"In a coordinate system with origin at the earth's center and axes 
always parallel to lines which are at rest with respect to the average 
positions of the fixed stars, that is, with respect to an inertial frame; 
see, for example, ref. 2, ch. I, or ref. ch. X. 
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gravitational field, and so forth. Let p denote the position vector 
of the target's center of mass. Then, 


P = 



( 2 ) 


where a is the position vector of the attack vehicle relative to the 
target . 

Using equation (2) to replace r* in equation (l) leads to 


d 2 a* 1 =* 


CM r» d 2 p* 


GM — > 
^ p 


1 — ^ 
- F 

m p 


(3) 


This differential equation is somewhat deceiving. The first two terms 
on the right side are "both large in absolute value, but for a « p (the 
most important case), the magnitude of their sum is comparatively small. 
Hence, for calculations, it is desirable to replace the entire sum on the 
right side of equation ( 3 ) by its approximate value derived as follows. 

Since, by equation (2), 

r = | p + a I 


it follows that 

r~3 = (p2 + 2p* . a* + a 2 ) (4) 

Consequently, for a « p, the approximate expression 

r~^ = p "^1 - 3 p" 2 P > ’ a 'j (5) 

is obtained by discarding terms of higher order in the series expansion 
of the right side of equation (4) . Another approximation which will be 
needed is the following. The last term on the right side of equation ( 3 ) 
represents the vector sum of the perturbing accelerations at the posi- 
tion r* of the attack vehicle and at time t. For a « p, it will be 
sufficiently accurate for present purposes to employ the approximation 


if = -i- F* ' 

m p m* p 


(6) 


where the primed quantities refer to the target (position p*). Since 
the exact equation of motion of the target's center of mass is 
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2_» 

d P _ CM -» 
dt 2 P 5 


i 

+ -==- F ' 
m' p 


(7) 


equations (5)^ (6), and (7) permit approximating the vector sum on the 
right side of equation (3) according to 


d 2 p* GM -» 
dt? " r3 


— F = 
m r P 


CM 


P 


(8) 


In view of the order of the approximation in equation (8), the third 
term on the left side of equation (3) is adequately approximated by 


M = Qi g* ( 9 ) 

r 5 p 5 

(Recall eq. (5)*) With the aid of the expressions (8) and (9), equa- 
tion (3) can be reduced to the approximate form 


d 2 a* 


dt 


-i? 

m p3\ 



( 10 ) 


Equation (10) Is, then, the vector equation of relative motion which 
applies accurately when the attack vehicle is near the target. Even in 
the case of zero thrust, the task of finding the general solution of 
equation (10) appears to be extremely difficult if for no other reason 
than that p cannot be expressed as a closed-form function of time 
when the target's orbit is noncircular . 


APPARENT GRAVITATIONAL ACCELERATION 


It is of Interest to consider the effect of the gravitational term 
in equation (10) for extreme cases: 


For case 1, 


— > - 

p . a = 0 


^ a* perpendicular to p*) 


In this case equation (10) can be written 

dt 2 m p? 


That is, the effect of gravity is to tend to reduce the distance a 
separating the attack vehicle and the target; moreover, this tendency 
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approaches zero as the separation a approaches zero. The analogy to 
the linear restoring force in the case of simple harmonic motion may 
also be noted. 


For cases 2, 

■ ) 

p . a = ±pa 


In these cases , 


same direction, same sense 1 
same direction, opposite sense] 


as 




and 




Equation (10) then reduces to 


d 2 i* 



I ?*+ 2 £§ 


That is , the effect of gravity is to tend to L increase the distance a 
separating the attack vehicle and the target; this tendency approaches 
zero as the separation a approaches zero . 


Let T denote the gravitational acceleration vector in equation (10). 
In the general case its magnitude 1b 



Consequently, 


GMa < _ < _CMa 


( 12 ) 


which shows the range of the magnitude of the gravitational acceleration 
in equation (10); moreover, it shows that case 1 and cases 2 are extremes. 
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The inequality on the right in expression (12) shows that when the attack 
vehicle is within a certain distance, say a Q , of the target, the magni- 
tude T of the gravitational acceleration term in the equation of relative 

motion (10) is less than 2 (t24 /p 2 ) a Q J P . This upper bound on T is a 
particularly useful check when a simplified a n alysis is performed with the 
effect of gravity neglected. 


The expression for the vector T shows that it lies in the plane of 

the two vectors a* and p*. (ihe special cases in which T and a* are 
collinear are case 1 and cases 2, discussed previously.) Having obtained 

the expression (ll) for the magnitude of it is desirable to have 

expressions for the angles 0 and 0 which r makes with a 

aT pr 

and p, respectively. Since 

Li. cos e 

pa ap 


where 


0 

ap 


is the angle which a makes with 


P> 


it follows that 


and 


cos 0 


aT 


a . r 3 cos2e ap - 1 
^ ||l + 3 cos 2 0 ap 


cos 0 

pr 



2 cos 9 ap 

l/l+5 cos 2 8 ap 


(15) 


m 


Equation (13) shows, for example, that T is perpendicular to a* for 

angles 0 ap = cos - that is, 54°1V and 125 °l6' • For both cases, 

the value of T is, by equation (ll), 1 / 2 (GM/p 2 )a/p. The corresponding 

angles 0 „ are, by equation (l4) , acute and obtuse , respectively. 

These cases, together with cases 1 and 2 described previously, are illus- 
trated schematically in figure 2. The acceleration field shown in fig- 
ure 2 may be thought of as being rigidly attached to the head of the 

changing vector pj the intensity of this attached field varies, of course, 
as l/p^, as well as varying directly as the distance a. Even for a tar- 
get in a circular orbit, in which case P = |f?| = Constant, the accelera- 
tion field shown in figure 2 rotates about the target in any coordinate 
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system whose origin is at the target and whose axes remain parallel to 
lines fixed in an inertial frame (cf. discussion that follows regarding 
proper representation of a*). The preceding theoretical discussion is 
amplified in part by figures 3 to 6. Figures 7 and 8 have been added 
to facilitate numerical estimates. The figures are self-explanatory. 

In the application of equation (10) it must be realized that any 
coordinate system to which the position vector a* may properly be 
referred must not only have its origin at the target's center of mass, 
but also have its axes always parallel to lines fixed in an inertial 
frame. This is necessary in order to be consistent with equations (l), 
(2), and ( 3 ) and the original definition of the vector r*. If r* had 
been referred to a coordinate system fixed in the earth, then terms 
corresponding to the centrifugal and Coriolis accelerations would have 
entered equation (l) because of the earth's rotation; the subsequent 
analysis would have been slightly more complicated. Equation (10) applies 
to the terminal phase of interception; consequently, the more complicated 
analysis is clearly undesirable when the tracking and guidance mechanisms 
for the terminal phase are totally contained in the attack vehicle . 

The apparent gravitational acceleration f* is mathematically identi- 
cal with a tide-generating acceleration and with the disturbing effect of 
a third body in celestial mechanics. This and some related matters are 
discussed in the appendix of this paper. 


TERMINAL GUIDANCE 


Guidance in the terminal phase of flight must be based at least upon 
the vector a or its variation with time, or both. A method which will 
be considered here Is somewhat similar to the familiar proportional 
navigation course (see fig. 9) snd is defined by the equation 


(o>0) (15) 

where c is the correction coefficient. Equation (15) states that the 
thrust F applied is to be perpendicular to a*, In the direction tending 
to reduce the rate at which a* is turning, and proportional to the 
angular rate at which a* Is turning. For this method of interception, 
equations (10) and ( 15 ) give 


d a c d /a*\ 

it 5 m 


+ 



= .0 


( 16 ) 
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as the equation of motion of the attack vehicle relative to the target 
in the terminal phase of flight. 


If, because of rocket operation, the mass m of the attack vehicle 
decreases appreciably during the terminal phase of flight, an additional 
equation must be solved simultaneously with equation (l6) . Let the exit 
speed of the propulsive exhaust gases, relative to the attack vehicle, 
be denoted by v ex ; then 


dm _ F 
dt “ v ex 


(IT) 


where F is the magnitude of the thrust vector F. The additional 
equation required follows from equations ( 15 ) and (l7)j that is, 


dm _ 
3t " 



( 18 ) 


This must be solved simultaneously with equation (l6) . If the correction 
coefficient c is a function of one or more variables, another equation 
expressing this dependence must be included. However, cases in which 
both m and c are practically constant are clearly of first importance 
and consideration will be restricted to these in the following discussion. 


SIMPLIFIED ANALYSIS NEGLECTING GRAVITY 


As in the case of equation (10), the gravitational term in equa- 
tion (l6) will generally make it necessary to resort to numerical methods 
of integration in precise calculations. Fortunately, however, the gravi- 
tational term approaches zero as the attack vehicle approaches the target; 
consequently, much can be learned by examining the simpler differential 
equation 


da c d i*_ 

dt 2 m dt a 


0 


( 19 ) 


which results upon setting f* equal to zero. Integrating equation (19) 
once gives 


da* , c a* _ / da*\ , c/a*\ 
dt m a \dt J. L m\a 


( 20 ) 


where the subscript 1 refers to values at the start of the terminal 
phase of flight. Since the vector on the right of equation (20) remains 
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constant in magnitude and direction, and since the vector 


c a* 


m a 


(c > 0) 


has the constant magnitude c/m, the attack vehicle is assured of hitting 
the target provided 


£ 

m 



( 21 ) 


which is easily seen from the vector diagram (fig. 10) corresponding to 
equation (20) . Let v^ denote the magnitude of the velocity vec- 
tor (da/dt ^ and let 0 be the angle which the latter makes with the 
negative of the vector a*. Then the inequality (21) is equivalent to 


or 




sin 0^ 


) 


2 


c > 1 V 1 
m 5 cos 0^_ 


(-f< 0 l<§) ( 22 ) 


( The inequality would be reversed for cos 0-^ < 0, corresponding to 
c < 0, but this does not result in a hit. ) This, then, is a sufficient 

condition for the attack vehicle to hit the target when gravity can be 
neglected and when the method of terminal-phase flight is that described 
analytically by equation ( 15 ) • 


The relative motion corresponding to equation ( 19 ), or the equivalent 
equation (20), is confined to the plane of the vectors (da/dt^ and 

(a/a^ or, in the trivial case for which these are collinear, the motion 

is confined to the line of the unit vector (a/a)^. Since, in accordance 

with equation ( 15 ), the thrust has zero component along a*, it follows 
from analytical mechanics (see, for example, ref. 2, ch. I, sec. 19 , or 
ref. 3, ch. V, sec. 1) that the radial component of the acceleration 
vanishes; that is, 




( 23 ) 


Since a ^ 0, equation ( 23 ) yields the relation 
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(24) 


The equality applies in the trivial case of straight-line motion. Assume 
now that relation (22) is satisfied, so that a hit will occur. Then, 
relation (24) shows that the magnitude -da/dt of the radial velocity 
is nonincreasing with time. This result is used in the following analysis 
of the variation of thrust magnitude during the terminal phase of flight. 


Hie identity 


d a? _ 1 /da* _ d& a*\ 
dt a a \dt dt a / 

and equation ( 19 ) yield 

d^ a* _ _l(c , o da\ d a* _ 1 d^a a* 
dt 2 a “a\m dt/cFt a ” a dt 2 a 


(25) 


The right side of equation (25) is the sum of two perpendic ul ar vectors; 
consequently. 


d d a* 
dt dt a 



ifc. 2 da\ |_d_ & 

a\m at/|dt a 


Therefore, by using equation (15) , 


dF 

dt 


2 

dt/ 


(26) 


is obtained for the time rate of change of the thrust magnitude F. 
Assume, as in the discussion following relation (24) , that relation (22) 
is satisfied. Hien, if the relation 


— = 2v n cos 8_ 
m 1 1 


( 27 ) 
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is also satisfied, equation (26) shows that the thrust magnitude F will 
not increase during the terminal phase of flight. 

Thus, if c/m satisfies the relations (22) and (27), the ma xlmim 
thrust required of the attack vehicle in the terminal phase occurs at 
the start, and its value F^_ is given by 


F 1 _ c V 1 l sin 9 ll 

m “ m a^ 


(28) 


If 0^. = 75° > relations (22) and (27) are satisfied by c/m = crv^ where 
a ^ 2, in which case equation (28) becomes 



( 29 ) 


NUMERICAL EXAMPLE FOR TYPE I INTERCEPTION 


Consider a type I vehicle for which the terminal phase of flight 
starts at a range of 50,000 to 100,000 feet with a relative speed of 
500 to 1,000 ft/sec and a heading error 0^ < 20°. Assume that o - 2 

(that is, c/m = 2v 1 ). The maximum required thrust F-^ given by equa- 
tion (29) is 

= 2 sin 20° = 13-7 ft/sec 2 

m 50000 

If the terminal phase of flight starts at time t^ and the target 
is hit at time t2, then (approximately) 


t 


2 



( 50 ) 


In the numerical example just considered, therefore, the maximum time 
possible in the terminal phase is approximately 




100000 
500 cos 20° 


= 213 sec 
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In this t ime interval a target moving at 25,000 ft/sec will have traveled 


(213) (25,000) = 5.33 X 10 6 ft = 1.01 X 10 5 miles 

This again is maximum for the numerical example. The corresponding max- 
imum fractional rate at which mass is expelled in propulsion of the 
attack vehicle is 


1 dm _ F l/ m 
m dt v ex 


15-7 

10000 


= 1.37 x 10" 5 


sec 


-1 


(31) 


where the exit speed v gx of the propellant gases has been taken as 

10,000 ft /sec. Thus, a substantial mass loss may occur in the extreme 
case. Finally, relations (12) show that the magnitude T of the 
apparent gravitational acceleration experienced by the attack vehicle 
during the terminal phase lies in the range 


0 < T < 2 ( 32 . 2 ) ft /sec 2 

- ■ 4000 


or 


0 £ r % 0.3 ft/sec 2 

Thus, the maximum value of T is fairly 'small in comparison with the 

maximum value F^/m = 13-7 ft/sec 2 and can be neglected in the present 
approximate calculation . 

Although the numerical example examined indicates that the applica- 
tion of an interception system of type I has attendant difficulties, 
considerable latitude exists for the adjustment of design parameters 
and none of the problems appears* to be extreme in nature . Consequently, 
this type of interception system is considered feasible. 


NUMERICAL EXAMPLE FOR TXPE II INTERCEPTION 


For interception of type II the relative speed v^ is of the order 

of 50,000 ft/sec. If a is taken as 2 (as before), equation (29) 
shows that reasonable values of Fj_/m imply large separation a-j_ and 

very small heading error 0^ at the start of the type II terminal phase. 
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For example, the values a-j_ = 500,000 ft and 9^ = 1° give 
F^/ni = 175 ft/sec^, T < 1.6 ft/sec ^ , and - t]_ = 10 sec (approxi- 
mately) . If it is assumed further that v ex = 10,000 ft/sec, the maxi- 
mum fractional rate at. which mass is expelled in propulsion becomes 
0.0175 sec"-'-. 

These numerical magnitudes serve mainly to indicate that the demands 
placed upon the type II interception system are severe. Even for a system 
which makes extensive use of ground equipment, the probability of a direct 
hit In this case appears small. A distance of nearest approach (miss 
distance) of the order of 1,000 feet will be assumed to be attainable in 
head-on interception. Furthermore, the timed dispersal of all or a part 
of the mass of the attack vehicle will be considered as a means of 
increasing the probability of a destructive hit to a value which is 
effectively unity. 


MASS DISPERSAL 


The dispersal of mass can be initiated by an explosion or some other 
means, which In the simplest case amounts to^ giving an impulse to each 
pellet or piece of matter ejected. 

For simplicity, assume that at time t the mass m is distributed 
uniformly over the surface of a sphere of radius 

b = V b (t - tj (32) 

where w^ is the radial speed of each piece of matter with respect to 
the center of mass (center of sphere) and t^ is the time at which dis- 
persal is Initiated. Let tg denote the time at which the target arrives 
at a point adjacent to the center of the sphere. At time t 2 the average 
surface density of mass on the sphere's surface is 

m _ m 

^b 2 2 4*[wb(t 2 - t i)] 2 

Hence, the masB m 6 which strikes the target is such that 


m B a 2A, 


m 

4*b 2 2 



B 
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d 


where Af is the target's frontal area. If m^ is the minimum 
destructive mass, then the probability of destruction -is a maximum when 


“tad ” 



(33) 


that is. 



Equation (33) has the following obvious and more important interpretation, 
however. Suppose the accuracy of the interception system is such that the 
center of mass of the attack vehicle is assuredly within a distance e of 
the target's flight path at time t 2 • Then, if s is substituted for b 2 

in equation (33) , the mi nlmimt mass m required to be dispersed is 

m = “md (3*0 


Thus, the dispersal speed w b and time interval t 2 - t^ are to be 
adjusted in accordance with 

b 2 = e = w b (t 2 - t ± ) 

If the attack vehicle Is approaching the target at a relative speed v^ 
at time t^ and if they are separated by the distance a^ at this 
instant, the distance at which dispersal is Initiated is approximately 

a i - (*2 - t i) v i - sy 

This distance will normally be held to a minimum so as to permit con- 
trolled flight of the attack vehicle for the longest possible time. 

The applicability of this type of interception is difficult to 
assess inasmuch as the minimum destructive mass depends upon the nature 
of the target, while the maximum error e depends upon the effective- 
ness of both the attack vehicle's control system and the ground equipment. 
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NUMERICAL EXAMPLE FOR MASS DISPERSAL 


In considering a numerical example it is convenient to rewrite equa- 
tion (34) in the form 

m _ 

m' m* kf 

where m* . is the target' s. mass . In view of the high relative speed, 
say v^ = 50,000 ft/sec, the ratio m^/m 1 will he taken as 10“5. 

Assume further that Aj = 100 ft^ and that e = 1,000 feet. Then, for 
these values of the parameters, 

JL. = io-5 6_x_io6 = 0<6 
m 10 2 

That is, the mass required to be dispersed is of the order of the mass 
of the target. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., July 23, 1958. 
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APPENDIX 

RELATION OF APPARENT GRAVITATIONAL ACCELERATION 
TO TIDE-GENERATING ACCELERATION AND TO 
DISTURBING EFFECT OF A THIRD BODY 


Hie apparent gravitational acceleration r experienced by the 
attack vehicle relative to a target moving in the earth's gravitational 
field is mathematically identical with the tide -generating acceleration 
field produced in the earth's oceans by the moon’s gravitational field. 
Superimposed upon the tidal acceleration field caused by the moon is a 
similar but somewhat less intense field due to the sun. Each of these 
tide -gene rating accelerations has the approximate form (cf. eq. (10)) 


p = _ GMo* 
p 5 

where M now refers to the mass of the tide -producing body (moon or sun), 
a* becomes the position vector of the ocean particle relative to the 
center of the earth, and p now refers to the position vector of the 
earth's center relative to the center of the tide-producing body. 

The expression (Al) for T will be derived from the approximate 
potential ( see , for example , ref . 4 or ref . 5 ) : 

u ■ I fs a2 (? ' oos2e «) . (A2) 

which Is well-known in tidal theory. For a given value of p, equa- 
tion (A2) defines a potential field about the point whose position vec- 
tor is p*. Thus, the independent variables In equation (A2) are a and 
0 a p, while p occurs simply as a parameter. The radial variable a and 
the colatitude 9 ap represent two of three possible spherical coordi- 
nates; the absence of the longitude angle cp in equation (A2) attests to 
the rotational symmetry of the potential field about the line of p*. 

Hence, in spherical coordinates, the gradient of the potential U is, 
by equation (A2) , 

ati 1 dU I du\ 

da* a d0 ap ; a sin 0 ap dq>/ 

= Tp(l - 5 cos 2 ©^, 3 sin 0 ap cos 0 ap , o) (A3) 

P 7 


grad U 


-■ 
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.The first or radial component of grad U in equation (A3) is in the 
direction of a* and can he expressed as 


GMa 

7 " 


(i 



(A4) 


The second component, which is in the direction of increasing 0 ap 

(see fig. 3), is perpendicular to a*, but can be expressed as the sum 
of a vector in the direction of a* and a vector in the direction of 
that is , 


(3 cos% - 

( 3 cos 6 

\ ap/a 

\ ap /P 


(A5) 


The third component of grad U in equation (A3) vanishes because of the 
rotational symmetry of the potential field. Adding expressions (A4) 
and (A5) gives 



Consequently, equations (Al) and (A6) show that 

T = -grad U 

which was to be proved. 

In celestial mechanics, the "disturbing effect of a third body" has 
a mathematical form which is identical with that of the tide-generating 
acceleration and the apparent gravitational acceleration of the present 
paper (see, for example, ref. 6, pp. 357-342). An example is the sun's 
disturbing effect on the moon's orbit about the earth. In this case r 
becomes the disturbing acceleration, M now denotes the mass of the 
sun, a* becomes the position vector of the moon's center relative to 
the earth's center, while p* becomes the position vector of the earth's 
center relative to the center of the sun. 


Thus, equation (Al) represents a first approximation for an accelera- 
tion in each of the three cases discussed, that is, terminal motion in 
Interception, tides, and lunar motion. In each application the ratio a/p 
is sufficiently small for the first approximation to be useful. 
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For a given value of the distance a, equation (Al) shows that the 
magnitude of the acceleration ? is mainly determined "by the value of 

the quantity Gm/p^, which will he termed the " disturbance parameter." 

In a given situation, such as the motion of a satellite about the moon, 
the disturbing effect of the earth can be compared with that of the sun 
by comparing values of the disturbance parameter for the two cases . This 
comparison shows that the disturbing effect of the earth on a satellite 
of the moon is nearly two hundred times as great as the disturbing effect 
of the sun on the satellite's motion relative to the moon. (This is in 
marked contrast to the GM/p 2 values for the two cases; the sun's 
attraction at the moon's position is more than twice the earth's. How- 
ever, each attracting body accelerates not only the moon's satellite 
but also the moon itselfj as a result, the disturbance parameter, which 
is a measure of the net acceleration of the satellite, varies inversely 
as the cube, rather than the square, of the separating distance p.) 
Therefore, successively better approximations for the motion of a satel- 
lite relative to the moon would include (in this order): (l) the moon's 

gravitational field, (2) the disturbing effect of the earth, and (3) the 
disturbing effect of the sun. 

Values of the disturbance parameter are given in table I to facili- 
tate the comparison of disturbing effects within each of five cases. In 
the first case, with which the present paper is concerned, the three 
values of the disturbance parameter show that the disturbing effects due 
to moon and sun are both negligible compared with that of the earth. The 
second (or earth- satellite) case is similar to the third (tidal) casej 
the values of the disturbance parameter here reflect the well-known fact 
that the moon's effect on the tides is somewhat greater than the sun’s. 
The fourth case has already been discussed, while the fifth is roughly 
analogous to the first case listed. 
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TABLE I.- VALUES OF THE DISTURBANCE PARAMETER 


© 

Disturbed 

body 

In motion 
near 

@ 

Disturbing 

body 

(mass M) 

© 

Approximate 
distance 
p between 

© and 
miles 

© 

Disturbance 

parameter, 

(24/p3, 

sec-2 

Attack 

Satellite of 

Earth 

>4 x 103 

<1.5 X 10 -6 

vehicle 

earth 

Moon 

*2.4 x 10 5 

*=8.6 x lO -14 



Sun 

«9.3 x 10? 

*=*3-9 x 10 " llf 

Satellite of 

Earth 

Moon 

2.4 x 10 5 

8.6 x 10 -11 *- 

earth 


Sun 

9.3 x 10 7 

3-9 x 10 - Ik - 

Ocean 

Earth 

Moon 

2.4 x io5 

8.6 x lo -11 *- 

particle 


Sun 

9-3 x 10 7 

3.9 x 10- lif 

Satellite of 

Moon 

Earth 

2.4 x 1C)5 

7-0 x 10" 12 

moon 


Sun 

9-3 x 10 7 

3-9 x lO -1 ^ 

Attack 

Satellite of 

Moon 

>1.08 x 10 5 > 

<9*5 x io-7 

vehicle 

moon 

Earth 

*2.4 x 105 

*»7-0 x 10 -12 



Sun 

*9-3 x 10 7 

*3.9 x 10- 114 - 



Point of initial detection and first 



Figure 1— Schematic diagram showing flight paths of two types of interceptor for the 

launching frcm the ground. 


case of 
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Numbers beside vectors indicate 
relative magnitudes 


54° 44' 





125° 16» 



54° 


Sphere, a ” Constant 


-Line of p* 


Earth's surface 


Figure 2.- Apparent gravitational acceleration T experienced by attack 
vehicle at small distance a frail target. (Acceleration field is 

rotationally symmetric about line of p* and varies in intensity as 


a 











Figure 5 .- Relation of angle 0^ which, f' makes with &) to angle 0 which a makes 

with p*\ (See eq.. (13).) 




9ap, dag 


Figure 6,- Relation of angle 


0 ( which 

pr v 

with p*). 


r makes with p ) to angle 
(See eg., (l4).) 



(which a* makes 




t. . [K a /. 

■‘max ” ^ ~ 5 r vxruni rexaGions (xk;)) 


p Q = Radius of earth - 5,959 miles 
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Q) Roc lost-applied acceleration , - £ ^ a, 

® in accordance Kith equation (l£) 

Direction of applied acceleration for 
a proportional navigation course 


Attack vehicle 



Figure 9«- Schematic comparison of direction of applied acceleration in 
present case (eq. ( 15 )} with that in the case of a proportional 
navigation course . 




